Experimental test of curvature-driven dynamics in the phase ordering of a two 
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We study electric field driven deracemization in an achiral liquid crystal through the formation 
and coarsening of chiral domains. It is proposed that deracemization in this system is a curvature- 
driven process. We test this prediction using the exact result for the distribution of hull-enclosed 
areas in two-dimensional coarsening in non-conserved scalar order parameter dynamics recently 
obtained [J.J. Arenzon et al., Phys. Rev. Lett. 98, 061116 (2007)]. The experimental data are in 
very good agreement with the theory. We thus demonstrate that deracemization in such bent-core 
liquid crystals belongs to the Allen- Cahn universality class, and that the exact formula, which gives 
us the statistics of domain sizes during coarsening, can also be used as a strict test for this dynamic 
universality class. 



Many aspects of the nonequilibrium relaxation of 
macroscopic systems still remain to be grasped. The do- 
main growth of two competing equilibrium phases after 
a quench from the disordered phase is a relatively simple 
out of equilibrium problem and in some cases the mecha- 
nism underlying coarsening is well understood (curvature 
driven, bulk diffusion, etc.) [l|. However, an important 
part of the description of these processes remains phe- 
nomenological and, to a certain extent, qualitative. 

According to the scaling hypothesis, a system in the 
late stages of coarsening is described by a scaling phe- 
nomenology in which there is a single characteristic 
length scale ("domain scale"), R(t), that grows with 
time. In consequence, all dynamical properties occur- 
ing on scales large compared to microscopic ones are de- 
scribed by scaling functions in which lengths are scaled 
by R(t). For example the pair correlation function, 
C(r,t) = (5(x,t)5(x + r,t)), where S takes the val- 
ues ±1 in the two equilibrium phases, has the scaling 
form C(r,t) = g[r/R(t)] [1]. Verifying the scaling hy- 
pothesis, and computing the scaling functions, has been 
a longstanding challenge. Recently, however, significant 
progress was made when the distribution rih(A, t) of hull 
enclosed areas (those enclosed by the outer boundaries of 
domains) was computed for scalar non-conserved order 
parameter dynamics in two dimensions and the scaling 
hypothesis verified for this quantity [2|, [3[ . The analyti- 
cally obtained distribution function was shown to be ro- 
bust and hold - to a high numericalprecision - in Monte 
Carlo simulations of the pure [2|, u| and disordered [4| 
bidimensional kinetic Ising Model (2dIM). The question 
remains as to whether more complicated experimental 
systems could also be described by such a universal for- 
mula. 

In this Letter we test this result experimentally in a 



liquid crystal system and we find a very good agreement 
with the theory. We thus demonstrate that the system 
belongs to the universality class of non-conserved scalar 
order parameter dynamics and that the exact formula is 
a universal property of these systems. In the following 
paragraphs we describe the experiment, and present a 
detailed analysis of the data. 

The experimental system chosen is one that exhibits 
electric field driven deracemization. Since the discovery 
of Louis Pasteur, more than 150 years ago, that chiral 
crystals can form from an achiral solution [5[, deracem- 
ization has been a fundamental question in the investi- 
gation of chirality. More recently, practical applications 
in, for example, drug design and synthesis, boosted re- 
search in this field, as the effect of most modern drugs is 
based on chiral molecules. Spontaneous deracemization 
in an achiral fluid system is very unusual and a topic of 
only recent interest [6|]. It can occasionally be observed 
in liquid crystalline systems [7|, |8|, [9[ , mainly formed by 
bent-core or so called "banana" molecules. The most 
likely reason for chiral conglomerate formation is steric 
interactions. This is also evidenced by computer simu- 
lations and theory, which indicate chiral conformations 
of on the average achiral molecules [l0[. Electric field 
induced switching between chiral domains was demon- 
strated in Ref. [Il[ . Kane et al [l2[ very recently exhib- 
ited the electric field driven deracemization of an achiral 
fluid liquid crystalline system, and gave a theoretical in- 
terpretation of the conglomerate formation in terms of 
a difference in the chemical potential of left and right 
handed molecules under electric field application. 

The liquid crystal employed in this investigation is 
comprised of a bent-core molecule, which together with 
cell preparation conditions is discussed in detail in 
Ref. [13[ . The studied cell has a gap of 5/im filled with the 
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liquid crystal, while lateral dimensions are much larger, 
approximately 1cm in each direction. We are thus effec- 
tively investigating a two-dimensional system. Domain 
coarsening was followed by temperature controlled po- 
larising microscopy (Nikon Optiphot-Pol microscope in 
combination with a Linkham TMS91 hot stage), with a 
control of relative temperatures to O.IK. Digital images 
were captured at a time resolution of Is with a pixel res- 
olution of N = 1280 x 960, corresponding to a sample 
size of 520 x 390/im 2 (JVC KY-F1030). Note that the 
imaging box is approximately 1/200 of the whole sample. 
Electric square- wave fields of amplitude E = 14V /im" 1 
and frequency / = 110Hz were applied by a TTi-TGlOlO 
function generator in combination with an in-house built 
linear high voltage amplifier. 

Cooling from the isotropic liquid, an optically isotropic 
fluid liquid crystal phase is formed, which exhibits no 
birefringence and thus appears dark between crossed po- 
larisers. The phase transition is first-order. On elec- 
tric field application, chiral deracemization occurs with 
domains of opposite handedness growing as a function 
of time. This coarsening process can easily be followed 
when the polarisers are slightly de-crossed by a few de- 
grees. The chiral domains of opposite handedness get 
larger as smaller domains disappear. At the same time 
the area distribution of domains with opposite handed- 
ness remains constant at an equal distribution of left- and 
right-handed domains, because there is an overall con- 
straint of zero chirality over the full sample that needs 
to be respected. This was checked experimentally and 
found to be true within experimental uncertainty due to 
finite size of the image window and thresholding. Still, 
such a global constraint (as opposed to a local one) is 
not expected to change the coarsening universality class 
which remains curvature-driven [l|. 

We performed 10 runs lasting lOmin each with pic- 
tures taken at intervals of 10 s on a single sample. Each 
run is initialized by heating the sample above the transi- 
tion temperature and subsequently cooling below it. The 
coarsening process in the low temperature phase is visu- 
alized in terms of domains, i.e., connected regions of the 
same handedness. In Fig. [T] we show a series of snapshots 
taken at times t = 0, 60, 120, 300 s. These pictures are 
then thresholded and an Ising spin Si is assigned to each 
pixel, where Si(t) = ±1 for pixels that belong to left 
or right handed domains, respectively. There are many 
spurious small domains that are related to the experi- 
mental system rather than to thermal fluctuations. The 
induced graininess is also reflected in the small r behav- 
ior of the pair correlation function C(r,t) and the small 
A behaviour of rih(A,t) as we shall see below. Still, at 
face value the domain geometry is the one of scalar non- 
conserved order parameter dynamics, as can be checked 
by comparing to the snapshots shown in [3[ for the 2dIM. 

The initial "magnetization density" in the imaging 
window of the liquid crystal, defined as the average of 
the spin variables over the box, m(0) = N' 1 5< (0), 
is not zero. This initial value is only approximately con- 




FIG. 1: The first snapshot displays the configuration right 
after the quench, t = s. The others are snapshots during 
the evolution t = 60, 120, 300 s. 



served by the dynamics, m(t) ~ m(0) = 0.2 ±0.1, but 
the actual value depends on the thresholding operation. 

We determine the growth law for the size, R(t), of typ- 
ical domains from a direct measure of the spatial corre- 
lation function, C(r,t) = jj Sill( s iCO s iCO )\\ri-rj\=r: 
The angular brackets indicates an average over the 10 
runs. The distance dependence of the pair-correlation 
at five equally spaced times, t = 100, . . . , 500s is dis- 
played with thin (red) lines in Fig. [2^a). As a conse- 
quence of the non-zero magnetization, C(r, t) does not 
decay to zero at large r. More strickingly, the curves 
are time- independent at distances r ~ 5 (C ~ 0.55) 
and they clearly depend on time at longer distances with 
a slower decay at longer times. Here and in what fol- 
lows we measure distances in units of the lattice spacing. 
The time-independence at short-distances and the long- 
distance decay are atypical, as can be seen by comparing 
to the spatial correlation in the 2<iIM displayed in the in- 
set to Fig. [2](a). We ascribe the lack of time-dependence 
at short scales and the further slow decay to the graini- 
ness of the experimental system. Indeed, in Fig. [2](a) we 
also show with thick dashed (black) lines the correlation 
in the 2<iIM where we have flipped, at each measuring 
instant, 10% of spins taken at random over the sample 
(the system dynamics are not perturbed and between 
measurements we use the original spins). By comparing 
the two sets of curves we see that the effect of the random 
spins is similar to the one introduced by the graininess 
of the system. This effect will also be important for the 
analysis of the hull-enclosed area distribution. 

The function C(r, t) obeys dynamical scaling, C(r, t) ~ 
g[r/R(t)]. We define the characteristic length-scale R(t) 
at time t by the condition C(R, t) =0.2 but other choices 
give equivalent results. The good quality of the scaling is 
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shown in Fig. [51(b). The time-dependence of the growing 
length R(t) is shown in the inset to Fig. [21(b) with points. 
The errorbars are estimated from the variance of the val- 
ues obtained from the 10 independent runs. We measure 
the growth exponent 1/z by fitting the long-time behav- 
ior of R(t), say fort > 30 s, and we find 1/z ~ 0.45±0.10. 
The exponent thus obtained is close to the theoretically 
expected value 1/2 for clean non-conserved order param- 
eter dynamics [1[ . The data suggest that for times longer 
than t ~ 30 s the system is well in the scaling regime. 



of hull-enclosed areas per unit area, rih(A, t) dA, with en- 
closed area in the interval [A, A+cL4], after a quench from 
high temperatures is 




FIG. 2: (Colour online.) Spatial correlation function at 
different times after the quench. (a) Experimental data 
at five equally spaced times, t — 100, . . . , 500s with thin 
(red) lines, and numerical simulation data in the 2dIM with 
10% randomly flipped spins at five equally spaced times, 
t = 100, . . . , 500 MCs with dashed (black) lines. We clearly 

notice the effect of graininess at very small scales (r ~ 5 in 
the experiment and r ~ 1 in the simulation), where there is 
no time-dependence in either case. Inset: the actual spatial 
correlation in the 2dIM. (b) Study of the scaling hypothesis, 
C(r,t) c± g[r/R(t)] in the liquid crystal, at the same times 
as in panel (a). Inset: the time-dependence of the growing- 
length scale. The slope of this line is 1/z ~ 0.45 =b 0.10. 

We now turn to the analysis of the distribution of hull 
enclosed areas. Each domain has one external perimeter 
which is called the hull The hull- enclosed area is the to- 
tal area contained within this perimeter. In [2, Q we de- 
rived an exact analytical expression for the hull enclosed 
area distribution of curvature driven two-dimensional 
coarsening with non-conserved order parameter. Using 
a continuum description in which the non-conserved or- 
der parameter is a scalar field we found that the number 



n h (A,t) = 2c h /(A + \ h tf 



(1) 



Ch = l/87r<\/3 is a universal constant that enters this 
expression through the influence of the initial condition 
and was computed by Cardy and Ziff in their study of 
the geometry of critical structures in equilibrium [14{ . 
Xh is a material dependent constant relating the local 
velocity v of an interface and its local curvature ft, in 
the Allen-Cahn equation, v = — (Xh/2ir)K |15l |. Equa- 
tion (pQ) can be recast in the scaling form rih(A,i) = 
(X h t)~ 2 f(A/X h t), with f{x) = 2c h /(x + l) 2 . In 
this way, scaling with the characteristic length scale, 
R(t) = V^ht, for coarsening dynamics with scalar non 
conserved order parameter in a pure system is demon- 
strated. The effects of a finite working temperature are 
fully encoded in the temperature dependence of while 
the same scaling function f(x) describes |3[ as sug- 
gested by the zero temperature fixed point scenario [l|. 




FIG. 3: (Colour online.) Scaling plot of the number density 
of hull-enclosed areas in: (a) the experiment; (b) the 2dlM 
with linear size L = 1280 and periodic boundary conditions 
evolving with non-conserved order parameter at T = 0. In 
the latter the measurements are done on a box with linear size 
£ = 1000. The lines are the prediction in Eq. (pQ). In the 2dIM 
case we exclude the spanning clusters from the statistics. In 
the insets we exclude all domains that touch the border while 
in the main panels we include them in the statistics. 

We counted the number of hull-enclosed areas in 
[A, A + gL4] to construct rih(A,t). It is important to note 
that the experimental data are taken using a finite imag- 
ing window and, in each sample and at each measuring 
time, many domains touch the boundaries (in contrast 
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to the numerical simulations in [2, 3, 4] in which we used 
periodic boundary conditions). The zero-chirality con- 
straint is not obeyed exactly, both because the image is 
just a subset of the whole sample and also because of the 
thresholding operation. 
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FIG. 4: (Colour online.) The number density of hull-enclosed 
areas for the liquid crystal sample at t = 100, 200, 400s, ex- 
cluding from the analysis the areas that touch the box border. 
Inset: nn for the 2d\M with 10% randomly flipped spins at 
measurement (see the text for an explanation) . In both cases 
the lines are the theoretical prediction (pQ) with ch = l/87r\/3 
and Xh = 2.1. 

In Fig. [31(a) we show the hull-enclosed area distribu- 
tion in the liquid crystal at three different times. In the 
main panel we included in the statistics the chopped ar- 
eas that touch the border of the image. The upward de- 
viation of the data with respect to the asymptotic power 
law A~ 2 is due to the finite image size. Indeed, domains 
that touch the border are actually larger but get chopped 
and contribute to bins of smaller A's and this induces a 
bias in the data. The inset displays nh removing from 
the statistics the areas that touch the border. The same 
anomaly appears in the 2dIM if one uses a finite imaging 
box within the bulk. To show this we simulated a system 
with L = 1280 and periodic boundary conditions and we 
measured in a finite square window with linear size 
i = 1000 using 100 independent samples. In Fig. 09(b) we 
show two sets of data for the 2dIM; in both cases we ex- 



clude the spanning cluster over the full system size. One 
set of data includes areas touching the border and lies 
above the theoretical curve. In the other set we elimi- 
nated these areas from the statistics and the datapoints 
fall on the analytic curve recovering the A~ 2 tail. 

The data in Fig. [H do not show any noticeable time- 
dependence at either small or large A. In the small area 
limit the time-independence can be traced back to the 
lack of time-dependence in the correlation function at 
distances r ~ 5 (which corresponds to A ~ irr 2 ~ 80), 
roughly the scale of the spatial graininess (see Fig. [2]). In 
the large area limit the time-dependence naturally disap- 
pears; structures with A ^> R 2 (t) are basically the ones 
already present in the initial condition and have not had 
time to evolve yet. In between these two limits the curves 
show a shoulder with a systematic time-dependence that 
is the most relevant part of our experimental data and 
it is very well described by the analytic prediction (pQ) 
shown with solid lines. To conclude we show that the 
random spins introduced by the measuring method are 
not only responsible for the time-independence of nh at 
small areas but also for the excess weight of the distri- 
bution in this region. In the inset to Fig. 0] we show 
the hull-enclosed area distribution in the 2dIM where we 
introduced 10% random spins at each measuring time. 
There is indeed a strong similarity with the experimen- 
tal data in the main panel that could even be improved 
by choosing to flip spins in a fine-tuned correlated man- 
ner. 

In summary, our experimental results for the hull- 
enclosed area distribution in the coarsening dynamics of 
the liquid crystal are in very good agreement with the ex- 
act analytic prediction for 2d non-conserved scalar order 
parameter dynamics presented in [2|, |3[ . 
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